Shell model description of normal parity 
bands in odd-mass heavy deformed 

nuclei 

C. Vargas^t J- G. Hirsch^t T. Beuschel^f J. P. Draayer^^ 

^ Departamento de Fi'sica, Centre de Investigacion y de Estudios Avanzados del IPN, 
Apartado Postal 14-740 Mexico 07000 DF, Mexico 
^ Instituto de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico, 
Apartado Postal 70-543 Mexico 04510 DF, Mexico 
^ Departament of Physics and Astronomy, Louisiana State University, 
Baton Rouge, LA 70803-4001, USA 

February 8, 2008 



Abstract 

The low-energy spectra and B(E2) electromagnetic transition strengths 
of i'^BEu, ^^^Th and ^^^^Dy are described using the pseudo SU(3) model. 
Normal parity bands are built as linear combinations of SU(3) states, 
which are the direct product of SU(3) proton and neutron states with 
pseudo spin zero (for even number of nucleons) and pseudo spin 1/2 
(for odd number of nucleons). Each of the many-particle states have 
a well-defined particle number and total angular momentum. The Hamil- 
tonian includes spherical Nilsson single-particle energies, the quadrupole- 
quadrupole and pairing interactions, as well as three rotor terms which 
are diagonal in the SU(3) basis. The pseudo SU(3) model is shown to be 
a powerful tool to describe odd-mass heavy deformed nuclei. 

PACS numbers: 21.60. Fw, 21.60.Cs, 27.70. -l-q 



The shell model is a fundamental theory that is applicable in nuclear, atomic 
and non-relativistic quark physics In its simplest formulation it provides a 
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natural explanation of magic numbers as shell closures and the energy spectra 
of closed shell ±1 odd-mass nuclei ^. Powerful computers and special algo- 
rithms for diagonalizing large matrices has allowed systematic studies of nuclei 
of the sd-shell Q and p/-shell up to A = 56 |5| . New methods for solving large 
scale shell-model problems in medium mass nuclei have also been developed 
. A shell- model description of heavy nuclei requires further assumptions that 
include a systematic and proper truncation of the model space 

In light deformed nuclei the dominance of the quadrupole-quadrupolc inter- 
action led to the introduction of the SU(3) shell model 0, and with it a very 
natural means to truncate large model spaces. Although realistic interactions 
mix different irreducible representations (irreps), the ground state wave func- 
tion of well-deformed light nuclei normally consists of only a few SU(3) irreps 
[H H, |l^, The strong spin-orbit interaction renders the usual SU(3) scheme 
useless in heavy nuclei, but at the same time pseudo-spin emerges as a good 
symmetry 

Pseudo-spin symmetry refers to the experimental fact that single-particle 
orbitals with j = 1 — 1/2 and j — — 2) + 1/2 in the shell 77 lie very close 
in energy and can therefore be labeled as pseudo spin doublets with quantum 
numbers j = j, fj = rj — 1 and 1 = 1 — 1. The origin of this symmetry has been 
traced back to the relativistic Dirac equation |l^, The pseudo SU(3) 
model capitalizes on the existence of pseudo-spin symmetry. 

In the simplest version of the pseudo SU(3) model, the intruder level with 
opposite parity in each major shell is removed from active consideration and 
pseudo-orbital and pseudo-spin angular momentum are assigned to the remain- 
ing single-particle states. The coupling of a deformed rigid-rotor core with one 
extra particle in a pseudo SU(3) orbital has been used to describe rotational 
bands and electromagnetic properties of heavy odd-mass nuclei and identi- 
cal normal and superdeformed bands fioll . 

A fully microscopic description of low-energy bands in even-even nuclei has 
been developed using the pseudo SU(3) model. The first applications used the 
pseudo SU(3) as a dynamical symmetry, with a single SU(3) irrep describing 
the whole yrast band up to backbending [^o| . A comparison of quantum rotor 
and microscopic SU(3) states provided a classification of the SU(3) irreps in 
terms of their transformation properties under tt rotations in the intrinsic frame 
and led to the construction of a operator which plays a crucial role in 
the description of the gamma band p3| . 

On the computational side, the development of a computer code to cal- 
culate reduced matrix elements of physical operators between different SU(3) 
irreps HJ] represented a breakthrough in the development of the pseudo SU(3) 
model. For example, with this code it is possible to include pairing, which is 
an SU(3) symmetry breaking interaction, in the Hamiltonian and exhibit its 



close relationship with triaxiality 1 25 , 2^. Full-space calculations in the p/-shell 
[ p7| in an SU(3) basis |ll| show that for a description of the low-energy spec- 
tra of deformed nuclei the Hilbert space can be truncated to leading irreps of 
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the quadrupole-quadrupole and spin-orbit (or pseudo spin-orbit) interactions. 
However, the inclusion of a pairing-type interaction is essential for a correct 
description of moments of inertia in such a truncated space. 

Once a basic understanding of this overall structure was achieved, a powerful 
shell-model theory for a description of normal parity states in heavy deformed 
nuclei emerged. For example, the low-energy spectra of many Gd and Dy iso- 
topes, their B(E2) and B(M1) transitions strengths for both their scissors and 
twist modes [ p8| and their fragmentation were successfully described with a 
realistic Hamiltonian |p9|| . 

In the present letter we introduce a refined version of the pseudo SU(3) 
formalism which uses a realistic Hamiltonian with single-particle energies plus 
quadrupole-quadrupole and monopole pairing inteactions with strengths taken 
from known systematics. The model is applied to three odd-mass rare earth 
nuclei: ^^^Eu, ^^^Tb and ^^^Dy. The results represent a full implementation of 
the very ambitious program implied in first applications of the pseudo SU(3) 
model to odd-mass nuclei performed nearly thirty years ago Jl^ . 

Many-particle states of Ua active nucleons in a given normal parity shell 7]a, 
a — v or TT, can be classified by the following chains of groups: 

{!""} {/a} {/a} 7a (Ao,Ma) ^ai^a 

U{n^) D U{n^/2) X U{2) D SU{3) x SU{2) D 

SO{-i) X SU{2) D SUj{2), (1) 

where above each group the quantum numbers that characterize its irreps are 
given and 7q and Ka are multiplicity labels of the indicated reductions. 

The most important configurations are those with highest spatial symmetry 
This implies that S-j^^i, = or 1/2, that is, only configurations with 
pseudo spin zero for even number of nucleons and 1/2 for odd number of nucleons 
are taken into account. 

We will describe ^^^Tb as a first example. It has 15 protons and 12 neutrons 
in the 50-82 and 82-126 shells, respectively. The number of nucleons in normal 
{N) and abnormal {A) parity orbitals is determined by filling the Nilsson levels 
with a pair of particles for /3 ~ 0.25 in order of increasing energy. This gives 

nf = 9, ni^Q, ^ 8, < = 4 (2) 

After decouphng the pseudospin in Eq. we get {/jr} — {2'*1}, {f^} = {2^} 
with = 1/2 and Si^ = 0. Table I lists the 15 pseudo SU(3) irreps, with the 
largest value of the Casimir operator C2, which were used in this calculation. 

Table I 
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The Hamiltonian contains spherical Nilsson single-particle terms for pro- 
tons and neutrons (-ffsp.7r[i/])j the quadrupole-quadrupole {Q ■ Q) and pairing 
{Hpair.T![v]) interactions as well as three 'rotor-like' terms which are diagonal in 
the SIJ(3) basis: 



C2. 

The term proportional to Kj breaks the SU(3) degeneracy of the different K 
bands |Q , the term represents a small correction to fine tune the moment of 
inertia, and the last C2 term is introduced to distinguish between SU(3) irreps 
with A and /i both even from the others with one or both odd ]2^ . 

The Nilsson single-particle energies as well as the pairing and quadrupole- 
quadrupole interaction strengths were taken from systematics 0; only a 
and b were used for fitting. Parameter values are listed in Table II and are 
consistent with those used in the description of neighboring even-even nuclei 

[!§. 

Table II 

Figure 1(a) shows the calculated and experimental js^ K — ^, | and ^ 
bands for ^^^Tb. The agreement between theory and experiment is in general 
excellent. The model predicts a continuation of the K = ^ band and over- 
emphasizes staggering in the K — ^ band. 

Figure 1 

The role played by each term in the Hamiltonian will be discussed in detail 
elsewhere In this letter we wish to emphasize that the pairing interaction 
is absolutely essential despite the strong truncation of the Hilbert space. To 
this end we present in part (b) of Fig. 1 the low-energy spectra of ^'^^Tb with 
the same Hamiltonian except that the pairing interaction has been turned off. 
It clearly exhibits the importance of the pairing interaction in building up the 
correct moment of inertia: the spectra without pairing is strongly compressed. 
It can also be seen that pairing affect the other energies in a similar way with 
an overall effect that resembles the introduction of a multiplicative factor in the 
Hamiltonian. We conclude that the proposed truncation scheme is justified and 
works as expected. 

Theoretical and experimental B(E2) transition strengths between yrast 
states in ^^^Tb are shown in Table III. The E2 transition operator that was used 



is given by |20 
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with effective charges Ctt = 2.3, Cu — 1-3. These values are very similar to those 
used in the pseudo SU(3) description of even-even nuclei |^ pof. They are 
larger than those used in standard calculations of B(E2) strengths ]30| due to the 
passive role assigned to nucleons in unique parity orbitals, whose contribution 
to the quadrupole moments is parametrized in this way. 

In Figure 2 (a) we present the low-lying energy spectra of ^^^Eu, including 
the K = I , I and ^ bands built with 7 protons in the normal parity subshell 
?7 = 3 and 8 neutrons in 77 = 4. There is a good agreement between the 
experimental and theoretical results. The model predicts a second state 
in the if = | band which is missing in the experimental spectra, as well as 
several other states in the excited bands. 

Figure 2 

It is interesting to notice that the ground state in ^^^Tb is |^ while in ^^^Eu 

it is I"*". Reproducing this effect is one of the successes of this theory; realistic 
single-particle-energies are required to get this ordering correct. 

The low energy spectra of ^^^Dy is presented in Fig. 2(b). There are three 
bands, with K — and ^ , respectively. As in the other cases the agreement 
between theory and experiment is remarkably good. In the K = ^ ground state 

band the ^ state is predicted to have an energy higher than the experimentally 
observed one. This departure of the experimental ground state band from the 
rigid rotor behavior may be related with a band crossing. The possibility of 
describing it by increasing the Hilbert space is under investigation. In the 
K = ^ band the | state lies higher than its | partner which contradicts the 
experimental results. As in the other cases, the model predicts several excited 
levels that are as yet undetected. 

It has been shown that normal parity bands in odd-mass heavy deformed 
nuclei can be described quantitatively using the pseudo SU(3) model. Only a 
few representations with largest C2 values and pseudo spin or 1/2 are needed. 
The Hamiltonian uses Nilsson single-particle energies, quadrupole-quadrupole 
and pairing interactions with strengths fixed by systematics, and three small 
rotor terms which with the others yield excellent results for energies and B(E2) 
values in A=159 nuclei. 

This work exhibits the usefulness of the pseudo SU(3) model as a shell model, 
one which can be used to describe deformed rare-earth and actinide isotopes by 
performing a symmetry dictated truncation of the Hilbert space. It opens up 
the possibility of a more detailed microscopic description of other properties of 
heavy deformed nuclei, both with even and odd protons and neutrons numbers, 
like g-factors. Ml transitions, and beta decays. 

This work was supported in part by Conacyt (Mexico) and the U.S. National 
Science Foundation. 
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Table Captions 



Table I: The 15 pseudo SU(3) irreps used in the description of "'^^^Tb bands. 
Table II: Parameters used in Hamiltonian ^. 

Table III: Theoretical and experimental B(E2) transition strengths for ^^^Tb. 



Figure Captions 



Fig. 1: Energy spectra of ^^^Tb, 'Exp' represents the experimental results 
and 'Theo' the calculated ones. Insert (a) shows the energies obtained with the 
Hamiltonian parameters listed in Table II, insert (b) shows the energies obtained 
without pairing. 

Fig. 2: (a) Energy spectra of ^^^Eu and (b) ^'^^Dy, with the same convention 
of Fig. 1. 
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Table 1: 



^^^Ell 0.00753 0l32 0)6 -0.0508 0009 0.0008 
159Tb 0.00753 0.132 0.106 0.0198 -0.0031 0.0008 
is^Dy 0.00753 0.132 0.106 0.0048 0.0006 0.0008 



Table 2: 
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Fig. 1 
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Fig. 2 
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